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Pair interactions whose Fourier transform is nonnegative and vanishes above a wave number Kq 
are shown to give rise to periodic and aperiodic infinite volume ground state configurations (GSCs) 
in any dimension d. A typical three dimensional example is an interaction of asymptotic form 
cos_K'or/r*. The result is obtained for densities p > pd where pi = Ko/2ii, p2 — {V^/8){Ko/tv)'^ 
and p3 — {1/8\/^){Ko/'!t)^ . At pa there is a unique periodic GSC which is the uniform chain, the 
triangular lattice and the bcc lattice for d = 1, 2, 3, respectively. For p > pd the GSC is non-unique 
and the degeneracy is continuous: Any periodic configuration of density p with all reciprocal lattice 
vectors not smaller than Kq, and any union of such configurations, is a GSC. The fee lattice is a 
GSC only for p > (l/6^/3)(7<'o/7r)^ 
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Crystallization of fluids is the paradigm of a continuous 
symmetry breaking. Its conceptual importance has long 
been recognized [1, 2], yet its derivation from first prin- 
ciples is still missing. Nature and computers can easily 
produce it, but a theoretical understanding of the emer- 
gence of a periodic order in continuous space as a result 
of a translation invariant interaction appears to be par- 
ticularly hard. The very first step along the way to a 
proof of a phase transition is to show that such inter- 
actions do have periodic ground states. This much more 
modest program has been advancing also very slowly, and 
for a long time the results were limited mainly to one di- 
mension [3-5]. A ground state configuration (GSC) is 
a minimizer, in a sense to be defined precisely, of the 
interaction energy. It is only recently that the mere exis- 
tence, without characterization, of infinite volume GSC 
was proved for a class of interactions in all dimensions 
[6]. The first two-dimensional example of an interaction 
giving rise to the triangular lattice as a GSC is even more 
recent [7] . The main concern of this paper is to provide 
examples of ground state ordering in three dimensions. 
The system we study is composed of identical classical 
particles interacting via pair interactions whose Fourier 
transform is nonnegative and decays to zero at a Kq < oo. 
Our results, although not predictive below a dimension- 
dependent density pd ~ Kq, are rather unexpected. At 
Pd a Bravais lattice (bcc for d = 3) is the unique periodic 
GSC. At higher densities the set of GSC is continuously 
degenerate: within certain limits, volume-preserving de- 
formations can be done on every GSC without cost of en- 
ergy, thus yielding other GSC. The degeneracy increases 
with the density in the sense that compressing any GSC 
results in a GSC of a higher density that can further be 
deformed. We can understand this proliferation of GSC 
as a consequence of the insensitivity of the interaction to 
details on a length scale shorter than K^^. That bcc lat- 
tice can be more stable than fee, should not surprise the 
reader. At equal densities the fee nearest neighbor dis- 
tance is slightly larger than the bcc one. For a purely re- 
pulsive interaction, the fee lattice is expected to be more 
stable; for a partly attractive interaction at some densi- 
ties the bcc lattice can have a lower energy. 



Definitions and notations. — We shall deal with trans- 
lation invariant symmetric pair interactions ip{r — r') = 
(p{r' — r). Rotation invariance will not be supposed. The 
A^-particle configurations {N < oo) are subsets of N 
points of M'* and will be denoted by Latin capitals B, R, 
X, Y. Infinite configurations with a bounded local den- 
sity will only be considered. The number of points in R 
will be denoted by 7V/j. If i? is a finite configuration, the 
interaction energy of R is U{R) = ^ J2r.r'eR,r^r' '^('"^ ~ 
r'). We wiU assume that and (^(k) = J^^^ (^(r)e"*''''' dr 
are absolutely integrable on R''. This ensures that both (p 
and (p{r) = (27r)~'* Jj^^ (,3(k)e''^ '" dk are continuous func- 
tions decaying at the infinity [8]. In the theorem below 

> implies U{R) > -ip{0)NR/2, hence, stability [9]. 
Let R and A" be a finite and an infinite configuration, 
respectively. The energy of R subject to the field created 
by X is given by 

C/(i?|X) = [/(i?) + ^^^(r-x). (1) 
reRxex 

Fix a real p. An infinite configuration X is called a grand 
canonical ground state configuration for chemical poten- 
tial p (pGSC) if it is stable against bounded perturba- 
tions, i.e., if for any bounded domain A and any R 

U{RnA\X\A)-pNRnA > U{XnA\X\A)-pNxnA (2) 

where X \ A is the set of points of X outside A. Be- 
cause every bounded domain is in a parallelepiped, these 
suffice to be considered. X is a canonical ground state 
configuration (GSC) if (2) holds for every R such that 
NrdA = AfxnA- Thus, any ^GSC is a GSC. If ip is su- 
perstable [9], the relation p i—f p (density) is invertible 
and any GSC is expected to be a pGSC for a suitable 
p. Local stability in the sense of Eq. (2) implies global 
stability, i.e., a GSC minimizes the energy density at the 
given density, cf. [10] and the end of this paper. 

A Bravais lattice B = {jyL=i^a^a\i^ & ^'^} will be 
regarded as an infinite configuration. Here are lin- 
early independent vectors and n = (ni, . . . , n^) is a d- 
dimensional integer. Any periodic configuration X can 
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be written as X ~ Uj^i{B + yj) where B is some Bra- 
vais lattice and B + y is B shifted by the vector y. For a 
given X, B is non- unique and we shall choose it so that 
J be minimum. Then, we call X a i3-periodic configu- 
ration. The reciprocal lattice is B* = {J2naha\n e Z''} 
where • b;3 = 2TTSaf3- The nearest neighbor distance 
(the length of the shortest nonzero vector) in B* will 
be denoted by qs*- This is related to the density via 
p{B) = CtypdqB* )'^, where ctypc is determined by the as- 
pect ratios and angles of the primitive cell of B. Let A 
be the parallelepiped spanned by the vectors La^a > ^ = 
{'^XaSia\0 < Xa < La}- We shall take La to be positive 
integers; then A is a period cell for B-periodic configura- 
tions, and the dual lattice A* = {^{ua/ La)ha.\n e Z**} 
contains B* . Next, we define the periodized pair interac- 
tion 



(3) 



and, for R in A, the periodized interaction energy 
UaIR) = 5 Lr.r'fEii.rT^r' '^A(r - r'). The sum defining 
<^A(r) is uniformly convergent, therefore <^a(i") is con- 
tinuous and tends to (p(r) as each La tends to infinity. 
Finally, let 
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where V{A) is the volume of A. /za tends to /x as A 

increases to W^. 

The main observation leading to the result presented 
below is as follows. Let ri, . . . , rjv be any finite configu- 
ration. If i^(k) > and is zero for |k| > Kq then 
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-r,) = {2nr' I ^(k)|Ve^'^-^pdk>0. 
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If we use periodic boundary conditions in a box contain- 
ing Vj , the integral is to be replaced by a sum. The k = 
term of this sum is structure-independent (for N fixed) 
and the rest is non-negative. Hence, any structure mak- 
ing the rest vanish is a GSC. But that is exactly what 
periodic structures accomplish, provided their shortest 
reciprocal lattice vector is outside the |k| = Kq sphere. 

THEOREM. Let both (p and (p he absolutely integrable 
in W^, <^(k) = (^(-k) > and (p{^^) = for |k| > Kq. 
We have the following results. 

(i) Let B he a Bravais lattice with qb* > Kq- Then every 
B-periodic configuration X is a GSC, and its energy per 
volume, e{X) = ^p[pip{0) — (p{0)], is minimum for the 
density p = p{X). On every period cell A, X n A min- 
imizes Ua{R) for fixed Nr = NxnA- Also, X creates a 
force-free field on test particles, i.e., U{r\X) is indepen- 
dent o/r. Any union of GSCs of the above type is a GSC 
(that can be aperiodic). 

(ii) There is a smallest density p^ at which qg- = Kq 
holds for a single Bravais lattice. If <^(k) > for < 



|k| < Kq, this B is the only periodic GSC. pi = KQ/2'!r, 
P2 = (V3/8)(Ko/7r)2 and = (1/8^2) (ifo/Tr)^, and 
the respective GSC are the uniform chain, the triangu- 
lar lattice and the hcc lattice. The fee, simple hexagonal 
(sh), simple cubic (sc) and hep lattices are GSCs at and 
above the respective densities pfcc = (l/6\/3)(ii'o/''')^, 
Psh = {V3/l6){Ko/7rf, psc = (l/8)(ii:o/7r)3 and phcp = 
(4/3V3)(i^o/7r)^ 

(Hi) Suppose, in addition, that (p{0) > and [p, + 
^Lp{0)]/(p{0) > pd. Then any B-periodic configuration 
X such that qs* > Kq (equivalently, p{B) > Pd) and 
p{X) = [p,+ |<^(0)]/(^(0) is a pbGSC, and its energy den- 
sity, e^(X) = e{X) — pp{X) = —^p{XY(p{Q), is mini- 
mum for the given p. If K is a period cell, AT n A mini- 
mizes Ua{R) — paNr. 

Because all the moments of (p are finite, (p is infinitely 
differentiable. A hard-core interaction can be added to (p 
provided that the close-packing density is larger than pd- 
Its only role is to restrict the set of allowed configurations. 
The set of GSC is reminiscent of a compressible fluid. 
The canonical and grand canonical ground state energy 
densities, 



^^[pm-m], 
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" 2"^'''"' "'^ 2^(0) 

are Legendre transforms of each other. 
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The p-dependence of Cp shows that the interaction is sta- 
ble if (p{0) > and is superstable [9] if <^(0) > 0. In the 
second case (7) expresses the cqiiivalcncc of the canon- 
ical and grand canonical ensembles. The corresponding 
density and chemical potential satisfy the equation 



p + <p{0)/2-(p{0)p = 0. 



(8) 



This linear relation breaks down for small densities be- 
cause p has to tend to —oo as p approaches zero. This 
implies a nonanalyticity, presumably at pd. 

Examples. — Based on a result on Fourier transforms 
[11], we can obtain fast-decaying (but not finite-range) in- 
teractions satisfying the conditions of the theorem. Take 
any locally integrable real function ^(k) = (?(— k) > 0, fix 
an£:>Oand define (p(k) = /n^,|^^^_^ 5(k')?7g(k— k') dk' 

where ?7e(k) = exp [—(1 — fc^/e^)~^] if fc < e and oth- 
erwise. This (f is infinitely differentiable. By inverse- 
Fourier transforming it we find (p to decay faster than 
algebraically. More interesting are the long-range inter- 
actions. They can be obtained by choosing (p> to be only 
finitely many times differentiable (namely, at |k| = i^o)- 
For example, in one dimension ipi^k) = Kq — \k\ for \k\ < 
Kq yields (p{x) = (1 — cosi4roa;)/7ra;^. In three dimensions 
rotation invariant examples can be obtained by starting 
with a function f{k) > such that / is three times con- 
tinuously differentiable and f{Ko) = f'{Ko) = 0. Then, 
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defining <^(k) = /(|k|) for |k| < Kq, by partial integra- 
tion 

= 2^ |[(fc/)"cosfcr]^'' - |^''"(fc/)"'cosfcrdfc| 

(9) 

For instance, f{k) = 7r^(fc + z){k + z){k — Kq)"^ with 
z = {Ko/10){l+3i) gives (p{r) = {13/lO)K'^cosKor/r^ + 
0{l/r^). The higher order terms make (p finite at the 
origin. One can verify that (p{0)/2(f{0) > p^, so this 
interaction has a continuous family of inequivalent ^GSC 
at /i = 0. Also, it has the bcc lattice as the unique 
periodic GSC at density (l/8V2)(A'o/7r)-\ 

LEMMA. Let A be a parallelepiped spanned by the vec- 
tors LaSicf Then 

Mr) = V{A)-' ^(k)e''^"'- (10) 

keA* 

In one dimension, for r = 0, Eq. (10) reduces to the 
Poisson summation formula [8]. 



Proof of the lemma. — The Fourier coefficients of ip^ 
are Jj^(pA{r)e~'^^''^ dr. Substituting the sum (3) for 
(pA{r), integrating by terms, and resumming, we obtain 
(^(k). The series on the right is a continuous periodic 
function whose Fourier coefficients are also f>(k). Be- 
cause of the completeness of the system {e*'' ''|k G A*} in 
the Banach space of integrable functions on A, equality 
for all r follows. 

Proof of the theorem. (i),(m) Consider the periodic 
configuration X = yjj^-^{B -\- yj) and let R be obtained 
from X by a bounded perturbation. Take a period par- 
allelepiped A of X large enough to contain the perturbed 
part, that is, R = X outside A. Recall that A* contains 
B* as a part. For a k in A*, 

J 

xexnA j=i 

where xb* (k) = 1 if k is in B* and otherwise. Using the 
lemma and Eq. (11), after some computation one finds 
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-piB) J2 ^(k)^e-^'^-'^^ Y ^ 

kGB* j=l rGfinA 



ik-r 



(12) 



Here p{B) = NBnA/V{A), the density of B. We can obtain U{X n A\X \ A) from Eq. (12) if we replace i? by X and 
reapply Eq. (11). Finally, the condition (2) reads 



xexnA 



ik-: 



E 

reflnA 



cik> (NRnA-NxnA) [m + ^(0)/2 - p(X)^(0)] 



+ p{B) Y ^(^) 



J J 
iVsnA l^e^''-"^ l^-^e-^'^''^ e^'^" • (13) 

j=i j=i reiinA 



Here we used iVxnA = JNbha, p{X) = Jp{B). If ((5(k) > 
and is zero for |k| > qs*, the left member is nonneg- 
ative, and the sum multiplying p{B) vanishes for all R. 
These and A^flnA = NxnA ensure that (13) holds true, so jf ^^^k) = for |k| > qe- then e{X) = \ (^(O)/?^ \Lp(fS)p 
that X is a GSC. If, moreover, p{X) = [^ + (p(0)/2]/(p(0), if also <^ > then X is a GSC and e(X) is the absolute 
the first term of the right member also vanishes for all minimum among configurations of density p. 
R, and X is a /iCSC. 

Choosing any period parallelepiped A, the canonical 
energy density of a B-periodic X of density p is Ua{R) = 



= \ p{Bf Y ^(k) I E e*'''' I ' - 



For ii in A and using again the lemma, 
1 



e(X) = 



1 



2ViA) 



2V{A) 



E E ^(---') 



(14) 
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+ 



Y <^a(x-x') - ^-^(0)^ 



2V^(A) 

Nr 
2V{A) 



E ^(k)lE^' 



ik-r I 2 



O^^keA* 
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x,x'eJcnA 



cf. Eq. (5). If i? = X n A, according to (11), the first 
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sum reduces to ^ k <E B* and vanishes completely if 
(^(k) = at |k| > qB' ■ When > 0, X n A minimizes 
Ua{R) for fixed Nr — NxnA and U\{R) — haNr, midcr 
the stronger conditions of (iii). One can also show the 
opposite implication: if X is periodic and X fl A„ mini- 
mizes Ua^ for a sequence A„ — > M** of period cells then 
X is a GSC. 

With A to be the primitive unit cell, 

t/(r|X)= E^(r-x)= ^A(r-x) (16) 

xex xexnA 

xexnAkeB* 

independent of r. This implies that the imion of GSCs is 
also a GSC that can be aperiodic; see [12] for details. 

(it) The one-dimensional case is obvious. In two di- 
mensions p{B) = (27r)~^|bi X b2|, in three dimensions 
p{B) = (27r)"^|(bi X b2) • h^]. We arc going to select B* 
by choosing in such a way that p{B) be minimum, 
on condition that qs* = min{| ^ naba|n ^ 0} = Kq. 
The resulting Bravais lattice at the given density will be 
denoted by B^. In two dimensions, let bi be one of the 
shortest vectors of B*, so that bi = qs*- Consider the 
shortest vectors of B* not collinear with bi . These form 
a star: with k, — k is also in the set. Choose b2 among 
them so that bi • b2 > 0. In the triangle of sides b\, 62 
and |b2 — bi| we have |b2 — bi| > 62 > so the largest 
angle is that of the vectors bi and b2. Therefore, this 
angle 0:12 is between 7r/3 and 7r/2. Under these restric- 
tions the minimum of |bi x b2| is obtained for 62 = hi 
and ai2 = 7r/3. These conditions define a triangular lat- 
tice for B2, and also a triangular lattice for B2, so that 
P2 = Ptr- lu three dimensions, the threshold densities of 
cubic lattices, at which the length of their shortest re- 
ciprocal lattice vector is Kq, can simply be computed. 
This gives the order pbcc < Picc < Pac and the values 
presented in the theorem. For the hexagonal lattice the 
threshold density depends on c/o, and the minimum is 
obtained for c/a = \/3/2. Its value, and that of the hep 
lattice at c/a = a/S/S and J = 2 are given in the theo- 
rem. It remains to convince oneself that the other Bravais 
lattices cannot give a smaller density. After some reflec- 



tion one can conclude that the conditional minimum of 
|(bi X b2) • bsl is attained by choosing fei = &2 = &3 = Kq 
and 7r/3 for the three angles between the three pairs of 
vectors. This specifies B^ as an fee lattice and -B3 as a 
bcc lattice. Thus, = pbcc- The last step of the proof 
of the theorem is to show that at the density pd no other 
periodic configuration X can be a GSC. Suppose there- 
fore that p{X) = Pd- If X is B-periodic but X ^ B then 
p{B) = pd/ J < Pd and therefore qB' < Kq. li X ^ B 
but B ^ Bd, then again qs' < Kq. Because i^(k) > 
for |k| < Kq, at least two nonzero vectors of B* give a 
positive contribution to the energy density (14). Hence, 
e{X) > ^pd[pd^{0)-^m = e{Bd) while p(X) = p{Bd). 
Then X is not a GSC because of the following. 

Let X and Y be two configurations of equal densi- 
ties but let e(y) < e{X). Then X cannot be a GSC 
of (p. In other words, every GSC minimizes the energy 
density at the given number density. Thus, no GSC 
can be metastable, cf. Ref. [10]. We only outline the 
proof. Take a large domain A of volume V such that 
NxnA = NynA- Estimate 

U{Xr\A\X\A) = U{XnA)+ J2 J2 ¥'(x-x'). 

xeJ>s:nAx'ex\A 

The first term is V{A)e{X) + o{V) while the second is 
o{V). Similarly, U{Y n A\X \ A) = V{A)e{Y) + o{V). 
Thus, U{Y n A|X \ A) < U{X n A\X \ A) for A large 
enough. This concludes the proof of the theorem. 

In summary, for a class of translation invariant pair 
interactions we have proved the existence of periodic 
and aperiodic ground states in classical particle systems. 
This result is probably the first of its kind in three di- 
mensions. Our finding, that ground states of innocent 
looking RKKY-type interactions such as cos Kor/r'^ are 
continuously deformable by volume preserving transfor- 
mations without cost of energy is quite unexpected. In 
another interpretation this means that existing crystal 
structures are stable against perturbations with interac- 
tions described in this paper. 
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